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1 4.  ABSTRACT  The  crust  of  southern  Eurasia  is  characterized  by  low  Q  values  (200-450)  for  1  -Hz  Lg  coda  everywhere  from  Spain  to  China.  The  north-south  extent  of  that  low-0  zone  varies 

fr°m  being  relatively  nareow  in  western  Europe  to  being  very  broad  from  the  Middle  East  to  China.  Those  variations  in  both  magnitude  and  width  first  became  apparent  with  the  development  of 
continent-wide  tomographic  maps  of  Lg  coda  Q.  The  most  recent  mappings  of  Lg  coda  Q  values  at  1  Hz  (&,),  completed  over  the  past  year,  show  that  the  four  regions  of  lowest  Q  values  occur  in 
regions  marked  by  high  levels  of  seismicity.  Qa  in  stable  regions  of  Eurasia  varies  between  about  450  and  900  with  the  highest  values  occurring  in  the  Indian  shield,  the  East  European  Craton  the 
Siberian  Craton  and  the  Khazakh  Platform.  The  maps  reveal  relatively,  and  unexpectedly,  low  values  of  Q0  in  the  Siberian  Trap  region  of  Siberia.  Variations  in  the  mapped  values  of  the  frequency 
dependence  of  Lg  coda  Q  at  1  Hz  (rj)  show  no  systematic  relationships  to  variations  of  Q0  in  either  southern  or  northern  Eurasia. 

Plots  of  point-spreading  functions — measures  of  how  well  features  can  be  resolved — reveal  that  features  with  dimensions  as  small  as  800  to  1000  km  can  be  resolved  The  standard  error  for  Q 
throughout  a  large  portion  of  Eurasia  is  50  or  smaller  and  between  50  and  100  in  most  other  regions.  An  exception  to  these  small  errors  occurs  in  a  region  that  extends  southwestward  from  Lake 
Baykal  where  they  are  as  large  as  200.  The  Baykal  values  probably  reflect  complexities  in  the  deep  crust  of  that  region.  Standard  errors  for  the  frequency  dependence  of  Q  also  show  anomalously  hich 
values  to  the  southwest  of  Lake  Baykal.  3  * 

Using  values  of  Qa  and  i\  obtained  for  Lg  coda  and  assuming  that  n  is  the  same  for  all  frequencies,  we  have  estimated  Rayleigh-wave  attenuation  coefficient  (  y*)  values  for  all  of  Eurasia  at  periods 
of  5,  10,  20  and  50  s.  In  our  studies  of  Rayleigh-wave  attenuation  and  its  regional  variation  for  various  periods  across  Eurasia,  we  have  assumed  that  the  values  of  Q0  and  rj  obtained  from  Lg  coda  can 
be  used  to  infer  values  of  shear-wave  Q  (QJ  and  its  frequency  dependence  (0-  In  that  process  we  assume  Q„  to  be  equivalent  to  an  average  value  for  shear-wave  Q  in  each  3°  x  3°  cell  of  Eurasia  and 
that  the  frequency  dependence  of  shear-wave  Q  can  be  obtained  using  an  empirically-derived  multiplicative  factor.  We  have  tested  this  assumption  for  several  crustal  shear-wave  Q  models  one  for  the 
Arabian  Peninsula,  one  for  the  Turkish/Iranian  Plateaus  and  seven  for  China  where  both  shear-wave  Q  models  and  Lg  coda  Q  information  are  available.  We  find  that  the  relations  Q  =  of  at  all 
depths  and  C  -  0.8 1\  for  the  depth  range  0-40  km  and  C=  0.5/j  for  depths  greater  than  40  km  provide  realistic  estimates  for  yR  in  the  period  range  5  -  70  s.  We  find  that  5-s  Rayleigh  waves°vary  between 
0  5  x  10  and  5.0  x  10  km' ,  10-s  waves  vary  between  0.2  and  2.2  x  10°  km  ',  20-s  waves  vary  between  0.1  and  1.1  x  10 km'1  and  50-s  waves  vary  between  0.06  and  0.3  x  lO  Mcm*1.  Sediments  play 
a  large  role  in  attenuating  5-s  and  10-s  Rayleigh  waves  in  some  regions  but  have  a  much  smaller  effect  for  20-s  and  especially  50-s  waves. 
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ABSTRACT 


The  crust  of  southern  Eurasia  is  characterized  by  low  Q  values  (200-450)  for  1-Hz  Lg  coda  everywhere  from  Spain 
to  China.  The  north-south  extent  of  that  low-£>  zone  varies  from  being  relatively  narrow  in  western  Europe  to  being 
very'  broad  from  the  Middle  East  to  China.  Those  variations  in  both  magnitude  and  width  first  became  apparent  with 
the  development  of  continent-wide  tomographic  maps  of  Lg  coda  Q.  The  most  recent  mappings  oiLg  coda  Q  values 
at  1  Hz  (QJ),  completed  over  the  past  year,  show  that  the  four  regions  of  lowest  Q0  values  occur  in  regions  marked 
by  high  levels  of  seismicity.  Q0  in  stable  regions  of  Eurasia  varies  between  about  450  and  900  with  the  highest 
values  occurring  in  the  Indian  shield,  the  East  European  Craton,  the  Siberian  Craton  and  the  Khazakh  Platform.  The 
maps  reveal  relatively,  and  unexpectedly,  low  values  of  Q0  in  the  Siberian  Trap  region  of  Siberia.  Variations  in  the 
mapped  values  of  the  frequency  dependence  of  Lg  coda  Q  at  1  Hz  (tj)  show  no  systematic  relationships  to  variations 
of  Q0  in  either  southern  or  northern  Eurasia. 

Plots  of  point-spreading  functions — measures  of  how  well  features  can  be  resolved — reveal  that  features  with 
dimensions  as  small  as  800  to  1000  km  can  be  resolved.  The  standard  error  for  Q„  throughout  a  large  portion  of 
Eurasia  is  50  or  smaller  and  between  50  and  100  in  most  other  regions.  An  exception  to  these  small  errors  occurs  in 
a  region  that  extends  southwestward  from  Lake  Baykal  where  they  are  as  large  as  200.  The  Baykal  values  probably 
reflect  complexities  in  the  deep  crust  of  that  region.  Standard  errors  for  the  frequency  dependence  of  Q  also  show 
anomalously  high  values  to  the  southwest  of  Lake  Baykal. 

Using  values  of  Qa  and  77  obtained  for  Lg  coda  and  assuming  that  77  is  the  same  for  all  frequencies,  we  have 
estimated  Rayleigh-wave  attenuation  coefficient  (ft)  values  for  all  of  Eurasia  at  periods  of  5,  10,  20  and  50  s.  In  our 
studies  of  Rayleigh-wave  attenuation  and  its  regional  variation  for  various  periods  across  Eurasia,  we  have  assumed 
that  the  values  of  Q„  and  l]  obtained  from  Lg  coda  can  be  used  to  infer  values  of  shear-wave  Q  (Q^)  and  its 
frequency  dependence  (£).  In  that  process  we  assume  Q„  to  be  equivalent  to  an  average  value  for  shear-wave  Q  in 
each  3°  x  3°  cell  of  Eurasia  and  that  the  frequency  dependence  of  shear-wave  Q  can  be  obtained  using  an 
empirically-derived  multiplicative  factor.  We  have  tested  this  assumption  for  several  crustal  shear-wave  Q  models, 
one  for  the  Arabian  Peninsula,  one  for  the  Turkish/Iranian  Plateaus  and  seven  for  China  where  both  shear- wave  Q 
models  and  Lg  coda  Q  information  are  available.  We  find  that  the  relations  Q„  =  Qac  at  all  depths  and  £=  0.8  77  for 
the  depth  range  0-40  km  and  £=  0.5  77  for  depths  greater  than  40  km  provide  realistic  estimates  for  ft  in  the  period 
range  5  -  70  s.  We  find  that  5-s  Rayleigh  waves  vary  between  0.5  x  10'3  and  5.0  x  10'3  km1,  10-s  waves  vary 
between  0.2  and  2.2  x  10'3  km'1,  20-s  waves  vary  between  0. 1  and  1 . 1  x  10'3  km'1  and  50-s  waves  vaiy  between  0.06 
and  0.3  x  10  km  Sediments  play  a  large  role  in  attenuating  5-s  and  10-s  Rayleigh  waves  in  some  regions  but  have 
a  much  smaller  effect  for  20-s  and  especially  50-s  waves. 
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OBJECTIVES 


Seismic  0  is  an  important  parameter  that  affects  how  rapidly  various  seismic  waves  attenuate  with  distance  as  they 
travel  through  the  Earth.  Its  value  can  vary  both  spatially  and  with  frequency.  The  overall  purpose  of  our  research  is 
to  better  understand  details  of  those  variations  throughout  southern  Asia  and  how  they  might  affect  our  ability  to 
detect  waves  generated  by  small  seismic  events,  to  estimate  magnitudes  and  yields  of  those  events,  and  to  determine 
if  they  are  explosions  or  earthquakes.  To  improve  understanding  in  those  areas  we  have  focused  on  tw  o  specific 
objectives.  Our  first  objective  has  been  to  expand  and  improve  previously  published  maps  of  regional  variations  of 
Lg  coda  Q  across  Eurasia  with  emphasis  on  regions  where  path  coverage  was  too  sparse  to  obtain  reliable  and  high- 
resolution  results.  Our  second  objective  has  been  to  use  those  new  maps  of  Lg  coda  Q  and  its  frequency  dependence 
at  1  Hz  ( Q0  and  rj ,  respectively)  to  develop  maps  of  Rayleigh- wave  attenuation  coefficients  for  various  periods 
(5-50  s)  throughout  southern  Asia. 

RESEARCH  ACCOMPLISHED 


Lg  Coda  Q  Variation 


We  presented  new  maps  of  Q0  and  77,  values  of  Lg  coda  Q  and  its  frequency  dependence  at  1  Hz,  for  all  of  Eurasia  in 
our  paper  for  last  year's  Seismic  Research  Review'  meeting  (Cong  and  Mitchell,  2005).  Both  last  year's  study  and 

this  one  use  the  method  developed  by  Xie  and 
Mitchell  (1990)  to  obtain  tomographic  maps  of  Q0 
and  77.  Details  about  our  methodology  can  be 
obtained  from  that  paper.  In  this  paper  we  again 
present  maps  of  Q0  and  77 .  This  repetition  of  the 
process  will  allow'  comparison  of  our  new  maps  of 
Rayleigh-wave  attenuation  variation  for  Eurasia 
with  the  Q0.  Also,  our  Q0  and  77  maps  have  been 
slightly  revised  because  of  a  bug  discovered  in  a 
new  code  that  we  used  to  determine  Lg  coda  Q  from 
individual  records. 


Figure  1.  Lg  coda  0at  1  Hz  (Q0)  for 


Figures  1  and  2  present  continent- wide  maps  of  Qc 
and  77  for  virtually  all  of  Eurasia.  They  were 
constructed  using  1383  determinations  of  those 
quantities  from  individual  seismic  records  for  paths 
confined  to  the  Eurasian  continent,  compared  with 


440  determinations  used  by  Mitchell  et  al.,  1997  in  an  earlier  study.  Data  coverage  was  best  in  regions  where 
seismicity  rates  are  highest;  i.e.,  in  the  southern  and  eastern  portions  of  the  continent.  The  coverage  was,  however, 
more  than  adequate  throughout  the  entire  continent. 


Figure  2.  Frequency  dependence  of 
Lg  coda  Q  at  1  Hz  (77)  for 
Eurasia. 


Figure  1  indicates  that,  with  the  exception  of  the 
Indian  platform,  Q0  is  generally  higher  in  northern 
than  in  southern  Eurasia.  But  even  in  northern 
Eurasia  unexpectedly  low  values  occur.  The  most 
conspicuous  of  these  is  in  central  Siberia  where  low' 
Q0  values  coincide  with  the  most  recent  mapping  of 
the  Siberian  Traps  (Reichow  et  al.,  2002).  The 
lowest  values  for  Q0  lie  in  the  four  regions  where 
seismicity  rates  are  highest.  These  are  the 
Kamchatka  Peninsula  in  northeastern  Asia,  the 
southeastern  portion  of  the  T ibetan  Plateau,  the 
Hindu  Kush  region  slightly  north  of  India,  and 
western  Turkey.  The  low  Q0  in  the  southeastern 
Tibetan  Plateau  extends  further  to  southeast  and 
covers  a  broad  region  of  enhanced  seismicity.  The 
Arabian  Peninsula,  even  though  it  is  a  stable 
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platform  over  its  entire  area  and  contains  a  shield  in  the  western  part  of  the  peninsula,  is  marked  by  00  values  that 
are  lower  than  most  other  stable  platforms. 

The  map  of  77  variation  in  Figure  2  shows  no  clear  relationship  to  0o  variation.  Low  77  correlates  with  low  Q0  in 
Kamchatka  in  northeastern  Siberia,  high  77  correlates  with  low  Q0  in  Spain  and  the  variation  of  the  large  low- 77 
region  in  the  north-central  part  of  the  map  appears  to  be  independent  of  Q0  variation  in  that  region.  This  differs  from 
some  other  continents  we  have  studied,  such  as  Africa  and  South  America,  where  high  Q0  values  are  associated  with 
low  77  values  and  low  Q0  values  with  low  77  values. 


Cong  and  Mitchell  (2005)  computed  point-spreading  functions  and  measures  of  resolution,  at  six  positions  in 
Eurasia.  They  showed  that  our  ability  to  resolve  features  in  the  crust  was  about  the 


»•  «*  o«  «a  11  It 

Relative  Vilar  of  rnf 


Figure  3.  Point  spreading  (unction  (psf)  plots  for 
six  selected  cells. 


same  everywhere  in  the  continent 
and  indicated  that  we  could  resolve  features  with 
dimensions  between  about  800  km  and  1000  km.  That 
paper  also  showed  that  the  standard  error  in  Q0 
determinations  was  less  than  50  in  most  places  and 
between  50  and  100  in  most  of  the  other  regions  of 
Eurasia.  A  region  to  the  southwest  of  Lake  Baikal 
exhibited  a  standard  error  greater  than  200.  Since  data 
coverage  in  that  region  was  good,  we  interpret  the  high 
standard  error  to  be  due  to  a  severe  lateral  variation  in 
crustal  properties  in  that  region.  Those  determinations 
appear  in  Figure  3.  Standard  errors  for  the  frequency 
dependence  of  Lg  coda  Q  at  1  Hz  (not  shown)  are  less  than 
0. 1  in  most  places  and  between  0. 1  and  0.2  in  most  of  the 
rest  of  Eurasia.  Like  Q0 ,  it  is  higher  to  the  southwest  of 
Lake  Baikal. 


Continent-Wide  Estimates  of  Rayleigh-Wave  Attenuation  Coefficients  (Yr)  for  Eurasia 


I 
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Figure  4.  Factor  by  which 
Qo  must  be 
multiplied  to  obtain 
a  good  empirical  fit 


Cong  and  Mitchell  (2005)  presented  estimates  of  Rayleigh-wave  attenuation 
coefficients  for  southern  Asia  for  periods  of  10,  20,  and  50  s.  In  this  paper  we 
extend  those  determinations  to  include  all  of  Eurasia  and  also  add 
determinations  for  all  of  Eurasia  at  a  period  of  5  s.  In  those  determinations, 
we  use  equations  developed  by  Anderson  et  al.  (1965)  to  estimate  Rayleigh- 
wave  attenuation  coefficients  from  Lg  coda  Q.  Lg  coda  Q  can  be  expressed 
by  the  equation 


where  Qa  is  Lg  coda  Q  at  1  Hz  and  rj  is  the  frequency  dependence  of  O  at  1 
Hz.  At  least  two  studies  (Herrmann  and  Kijko,  1983;  Campillo  et  al.,  1985) 
have  found  that  Q0  in  the  above  equation  is  a  good  approximation  of  the 
average  value  of  shear-wave  Q  ( Q in  the  crust.  An  expression  for  shear- 
wave  is 


Q.=Qf 


(2) 


where  £is  the  frequency  dependence  for  Qu.  We  take  Q^m.  this  expression  to  approximate  Q0.  We  had  initially 
thought  we  could  let  £ approximate  77  but  found  that  we  could  not  duplicate  measured  values  of  Rayleigh-wave 
attenuation  by  doing  that.  We  then,  by  trial  and  error,  empirically  obtained  a  depth-dependent  factor  that,  when 
multiplied  by  77,  provided  a  good  fit  between  ft  predicted  by  Lg  coda  Q  and  observed  ft  curves.  That  factor  is  0.8 
for  the  upper  40  km  in  the  Earth  and  0.5  at  greater  depths. 
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Equation  3  shows  the  equation  of  Anderson  et  al.  (1965)  as 
modified  by  Mitchell  and  Xie  ( 1 993)  to  take  into 
consideration  depth-variable  frequency  dependence  of  Q. 


r.= 


n 

cifr. 


v 


P, 


to. 


(3) 


Figure  5.  Kata  period  of  5  s. 


Figure  6.  Y«  at  a  period  of  10s. 


Y  il»‘ 


In  this  expression,  CR  refers  to  Rayleigh-wave  phase 
velocity,  T  to  period,  to  shear-wave  velocity  in  layer  /, 
to  compressional-wave  velocity  in  layer  1  and  to  shear- 
wave  Q  in  layer  /.  The  subscripts  to  the  right  of  the  partial 
derivatives  refer  to  quantities  held  constant  during  the 
differentiation,  p  in  those  subscripts  refers  to  density.  The 
Rayleigh-wave  phase  velocities  and  the  partial  derivatives  in 
equation  3  were  computed  for  every  3°  x  3°  cell  in  our  Q0 
map  using  many  shear-wave  velocity  models  for  Eurasia 
provided  by  M.  Ritzwoller  and  A.  Levshin. 

We  have  tested  the  utility  of  our  conversion  for  several 
crustal  shear-wave  Q  models,  one  for  the  Arabian  Peninsula, 
one  for  the  Turkish/Iranian  Plateaus,  and  seven  for  China 
where  both  shear-wave  Q  models  and  Lg  coda  Q 
information  are  available.  As  indicated  above  we  find  that 
the  relations  QR  =  Q0C  at  all  depths  and  0.8 tj  for  the 
depth  range  0-40  km  and  0.5  77  for  depths  greater  than  40 
km  provide  realistic  estimates  for  yR  in  the  period  range 
5-70  s.  We  find  that  5-s  Rayleigh  waves  vary  between  1.5  x 
10'3  and  5.0  x  lO^km*1,  10-s  weaves  vary  between  0.2  and 
2.2  x  10'3  km*1,  20-s  waves  vary  between  0. 1  and  1.1  x  10*3 
km'^and  50-s  waves  vary  between  0.06  and  0.3  x  10‘3  km'1. 
Sediments  play  a  large  role  in  attenuating  5-s  and  10-s 
Rayleigh  waves  but  are  less  important  for  20-s  and, 
especially,  50-s  waves. 


Figure  7.  Yr  at  a  period  of  20  s. 
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Figure  8.  %  at  a  period  of  50  s. 


At  both  5-and  10-s  periods  there  are  several  regions  where 
attenuation  is  quite  high.  Some  of  these,  such  as  the  Barents 
north  of  the  continent  and  the  Black  Sea  region  are 
obviously  associated  with  thick  accumulations  of  sediments. 
There  are  also  several  smaller  regions  with  high  attenuation 
rates.  We  checked  the  velocity  models  for  several  of  these 
regions  and  found  that  all  of  them  were  also  marked  by  thick 
accumulations  of  low- velocity  sediments.  As  periods 
increase,  the  smaller  regions  with  high  attenuation  disappear 
and  the  Q  variations  are  broader  with  fewer  small-scale 
features.  This  is  due  to  the  smaller  sensitivity  of  the  longer 
wavelength  waves  to  detail  and  to  the  fact  that  those  waves 
are  less  sensitive  to  near-surface  features. 

At  50-s  periods,  the  Rayleigh- wave  attenuation  coefficient 
map  resembles  the  1-Hz  Lg  coda  0  map.  This  is  probably 
because  the  Lg  coda  samples  the  entire  crust  and  is  not,  like 
5-  and  10-s  Rayleigh  waves,  dominated  by  the  upper  crust. 
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Similar  to  the  0  map  tor  Lg  coda,  the  50-s  Rayleigh-wave  attenuation  map  shows  regions  of  low-attenuation  (or 
high  0  in  the  three  stable  regions  of  northern  Eurasia,  as  well  as  the  Indian  platform.  It  also  shows  low  attenuation 
for  most  of  the  Arabian  Peninsula  and  southeastern-most  Asia  in  contrast  to  the  relatively  high  attenuation  (or  low 
0  exhibited  by  the  0  map.  This  occurs  because  the  frequency  dependence  parameters  (rj)  in  Figure  2  are  relatively 
low  in  both  of  those  regions. 


CONCLUSIONS  AND  RECOMMENDATIONS 

New  maps  ot  Lg  coda  Q  at  1  Hz  (0  and  rj)  have  been  used,  along  with  regionalized  seismic  velocity'  models,  to 
develop  continent-wide  maps  of  Rayleigh-wave  attenuation  coefficients  at  periods  of  5,  10,  20,  and  50  s.  The  Lg 
coda  Q  maps  were  found  to  resolve  features  between  about  800  km  and  1000  km  in  dimension  and  to  have  relatively 
low  standard  deviations  everywhere  but  a  region  southwest  of  Lake  Baikal.  The  Rayleigh-wave  attenuation 
coefficient  maps  were  determined,  based  upon  published  work,  assuming  that  0  in  each  3°x  3°  cell  of  Eurasia 
represents  an  average  value  for  the  shear-wave  Q  (0)  there  and  that  the  frequency  dependence  of  0  can  be 
obtained  from  an  empirically  determined  depth-dependent  multiplicative  factor.  That  factor  is  0.8  alt  depths  between 
0  and  40  km  and  0.5  at  greater  depths.  It  was  tested  for  paths  in  China,  the  Arabian  Peninsula,  and  the  Iran/Turkey 
Plateaus  and  was  found  to  be  satisfactory  in  all  cases  but  those  where  the  0  or  rj  may  not  be  well  determined.  We 
find  that  5-s  Rayleigh  waves  vary  across  Eurasia  from  0.5  x  10°  and  5.0  x  10’3  km'VlO-s  waves  vary  from  about  0.2 
tp  2.2  x  10  km  ,  20-s  waves  vary  from  about  0.1  to  1.1  x  10’3  km’1, and  50-s  waves  vary  from  about  0.06  to  0.3  x 
10  km’1.  Sediments  play  a  large  role  in  attenuating  5-s  and  10-s  Rayleigh  waves  in  some  regions  but  have  a  much 
smaller  effect  for  20-s  and  especially  50-s  waves. 

The  validity  of  our  empirical  relation  for  obtaining  shear-wave  Q  frequency  dependence  used  in  computing 
Rayleigh-wave  attenuation  coefficients  is,  so  far,  based  only  upon  Rayleigh-wave  attenuation  information  in 
tectonically  active  regions.  For  that  reason,  we  recommend  obtaining  Rayleigh-wave  attenuation  measurements  for 
regions  in  northern  Eurasia  that  are  relatively  stable.  Also,  our  determinations  assume  that  the  ratio  of 
compressional-wave  Q  to  shear-wave  Q  is  equal  to  2  in  crystalline  rock  beneath  sediments.  Recent  work  suggests 
that  that  ratio  is  closer  to  1,  at  least  in  tectonically  active  regions.  We  recommend  conducting  combined 
determinations  of  compressional-wave  Q  and  shear- wave  Q  in  a  variety  of  tectonic  regimes  to  obtain  values  for  that 
ratio  in  the  crust. 
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